Random walk on p— adics in glassy systems 
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We show that p— adic analysis provides a quite natural basis for the description of relaxation in 
hierarchical glassy systems. For our purposes, we specify the Markov stochastic process considered 
by S. Albeverio and W. Karwowski. As a result we have obtained a random walk on p— adic integer 
numbers, which provide the generalization of Cayley tree proposed by Ogielski and Stein. The 
temperature-dependent power-law decay and the Kohlrausch law are derived. 
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C/2 ; I. INTRODUCTION 

The growing interest in spin glassy dynamics has stimulated intense research in experimental [1-4] as well as 
theoretical domain [5-17]. The theoretical studies have been successfully developed in the following two directions: 

(i) scaling theory for growing domains and droplets [5-7] 

(—| , (ii) the hierarchical structures in ultrametric spaces for relaxation dynamics [8-17]. 

o 

q ' The goal of this paper is to investigate the second item (ii), using p— adic analysis. The temperature cycle experiments 
[1-4] are usually interpreted by the existence of a continuous hierarchy in spin glass dynamics. The hierarchical 
organization of metastable states depends on the temperature. It is assumed that at given temperature T the system 
occupies restricted part of the phase space inside one of the free energy valleys. After cooling down to temperature 
T — AT, where AT is rather small in the range 1°K — 2°K, such a part of phase space is divided into several smaller 
dcsccndcnt valleys separated by energy barriers. If the temperature is lowered the system has its new landscape of 
, valleys but the information about the previous states is not lost (memory effect). So, after heating the system to 
temperature T all descend states merge together into much smaller numbers of their ancestor states. The system 
in temperature T continues the relaxation process. After heating the system to critical temperature T c all occupied 
states merge into unique common ancestor state and phase space looks as in paramagnetic phase. 

From temperature cycling experiments it can also be seen, that the relaxation process at temperature T during time 
t is practically the same as the process at temperature T" = T — AT during some shorter time t' < t. The free energy 
barriers separating valleys at temperature T' = T — AT are higher than those at the temperature T and corresponding 
amount of time t needed to occupy the respective part of the phase space at temperature T is smaller than in the 
case of temperature T' . Hierarchical structure of metastable states corresponds to such structure expressed in terms 
of pure states in the sense of Parisi solution [17]. If we lower the temperature the energy barriers separating states 
O ■ are becoming higher, even infinite. 

The dominant process in the relaxation dynamics has been described as the hopping between the physical states 
[4]. Such mechanism leads to the studies of the models with a hierarchy-constrained dynamics with ultrametric 
topology. These models did appear very useful for the description of complex systems, called glassy systems, with 
highly degenerate metastable states. 

In this paper we consider dynamics for glassy systems using a random walk on p-adic numbers. We specify, to our 
purposes, the Markov stochastic process, considered by S. Albeverio and W. Karwowski [18, 19], for the description 
of dynamics in ultrametric spaces. We use the space of p— adic integers, which for our purposes is an appropriate 
example of ultrametric space. We should notice that the p— adic integers may be represented by the bottom of the 
regular, infinite Cayley tree. Let us recall that in Cayley tree of order p every branch at each level splits into p other 
branches. For the illustration we present on Fig. 1 the Cayley tree with p = 2 from which we obtain, after infinity 
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number of splitting, the 2-adic integers as describing the bottom of the tree. 
We can introduce a random walk in p— adic integers Z p analogously as a random walk in a ultrametric space 




FIG. 1: Cayley tree with p=2. The bottom represents Z2. 



represented by leaves of regular, finite Cayley tree. The pure physical states of the system are represented by these 
leaves of infinite Cayley tree. In order to construct the Markov process on p-adic integers we firstly consider Markov 
process on p-adic balls with the same finite radius and finally we contract the ball radius to zero. In that way we 
obtain transition probabilities in the space of p— adic integers which correspond to the ones considered by Ogiclski 
and Stein [9], in a case of finite ultrametric space. 

The plan of the paper is the following. In Section 2 we introduce basic mathematical notions and the properties 
of p— adic integers. We describe as well a Markov process on p-adic integers. In Section 3 we pass to the physical 
application. We consider the thermally activated random walk on ultrametric space described as a Markov process on 
the space of p-adic balls. The temperature-dependent power-law decay and the Kohlrausch law are derived. Let us 
notice finally that recently one observes an increasing interest in the application of p-adic numbers in mathematical 
physics [20]. The p-adic analysis was used to study stochastic processes [20, 21], especially to ultrametric jump 
diffusion [22-24]. 

II. CONSTRUCTION OF MARKOV PROCESS ON THE p-ADIC INTEGERS. 

Before we consider a Markov process with p— adic integers as a state space let us first introduce some notation and 
basic properties of p— adic numbers. More details may be fined, for example in [25] or [20]. 

A. p-adic integers 

Let p be an arbitrary prime number and let Z p denote the set of p— adic integers. A p— adic integer is a formal 
series J2i>o a «P* coefficients ai satisfying < <Xj < p — 1. With this definition, a p— adic integer a = X)i>o a iP l 
can be identified with the sequence (cij)j>o of its coefficients. 

In order to introduce a distance between p— adic integers a and b let us first consider an order of a p— adic integer. 
The order of a p— adic integers a = (dj)j>o is the smallest m for which a m ^ 



ord p (...a2 ai do) = min{« : dj 7^ 0}, 



(1) 



with the convention that maximum of the empty set is equal infinity. 
Notice that if x, y are p— adic integers then 

ord p (xy) — ord p (x) + ord p {y), (2a) 
ord p (x + y) > min(ord p (x), ord p (y)). (2b) 

Now, in term of the function ordp(-), we may introduce in the space of p— adic integers the norm || • || p 

\\x\\ p =p- ord *W, (3) 

and the p— adic metrics 

d p (x,y) = \\x-y\\ p =p-° rd ^-yl (4) 

It is easy to check that the function d p is a metric, and moreover due to inequality (2b), fulfils the following ultamctric 
condition 

d P (x,y) <max(d p (x,z), d p (z,y)). (5) 
For M > and a e Z p we may define a closed p— adic ball B(a, M) with center a and radius p~ M 

B(a,M) = {x e Z p : d p (x, a) <p~ M } . (6) 
If centre of the ball a has p— adic representation a = J2"jLa a jP ^ then the ball B(a, M) is completely determined by 

{o]m = «m, clm-i, . . . , a . (7) 
Let us recall some basic properties of the p— adic balls: 

(i) for any p— adic balls B(a, M) and B(b, M) we have 

B(a, M) n B(b, M) = or B(a, M) = B(b, M), (8) 
which means that any p— adic balls are disjoint, or one is enclosed in another. 

(ii) every point of the ball is the centre of this ball 

if a; € B(a, M) then B(x, M) = B(a, M). (9) 

(iii) each ball B(a, M) of radius p~ M may be represented as a finite union of disjoint balls S(a,, M + 1) of radius 



p -(M+l) 



p-1 

B(a,M)=\jB( ai ,M + l), (10) 

(iv) (iii) implies that for disjoint balls S(a,, M) of radius p~ M 

P M -i 

Z p = |J S(ai,M). (11) 

i=0 

Let us finally define the distance between two arbitrary balls B\ and B 2 as 

d p (B u B 2 ) = inf{d p {x, y) : x £ B u ye B 2 }. (12) 

With such definition two balls with representations {cl\m = a M, a M-i, ■ ■ ■ , a k, a-k-i, ■ ■ ■ , °o an d {b} M = 
&m, &m-i, • • • ,bk, ctfe-i; • • • , a o, differ for the first time at k — th position, we have 

d P ({a} M , {b} M ) =p- h - (13) 
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B. Random walk on p— adic balls 

Let {X(t), t > 0} be a Markov process with the state space Z p and transition rates between states depending only 
on p— adic distance between these states. To define such a process we a follow general construction given in [19]. In 
order to do this we represent Z p as a finite union of disconnected balls Bf with the radius p~ M 

P M -i 

z P = (J B ™ ( 14 ) 

i=0 

For a fixed M we may consider a finite state space Markov process {X M (t), t > 0} with the state space E M = 
{Bf, Bf , . . . , Bf Ml \, set of p— adic balls with the radius p~ M , and transition probabilities from ball Bf at time 

to ball Bf at time t, P^\t), defined as 

P^\t) = P(X M (t) = Bf\X M (0) = Bf) 



= P (X(t) e Bf I X(0) g Sf) . (15) 
Transition probabilities Pj;f\t), = 0, 1, . . . ,p M — 1) are solutions of the system of Kolmogorov equations 



d_ 

It 



0<k<p M -1 



with the initial condition -Pij (O) = <5y, where q\f^ the infinitesimal transition probability and q\ M ^ the intensity of 
stay in state i of the process are defined for any pair i, j of different states represented by balls Bf , Bj^ with radius 

p- M ,by 

(M) r t,j \ n ) 

q\ = hm—±- , (17 

and for any state i, by 

(M) ,. * _ ^ W /, N 

« = te — h — • (18) 

Observe, that each ball Bf is a union of disjoint balls Bf +1 , (k = 1, . . . ,p) of radius p~( M+1 ). We stress here, 
that depend only on p— adic distance between balls Bf and Bf . For this reason we may represent transition 
probabilities of the process Xm(<) by appropriate transition probabilities of the process Ijw+i(t) 

Plf\t) = P(X M (t) = Bf\X M (0) = Bf) 



P(^(t)e[jBf+fX(0)e[jBf k ^ 



\ fe=l k=l / 

= P P (X M+1 (t) = Bf +1 | X M+1 (0) = Bf +1 ) 

= pPg +1 \t)) (19) 

This equality leads to the recurrence relation for the local characteristics and qlf +1 ^ of the processes {XM(t)} 
and {X M +i(t)}. Let dist p (Bf ', Bj^) = p"", then we define 

gg° =p- M u(-M,M-n). (20) 
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we may notice that definition of q\ • ' and (19) impies 

u(-M + l,m-l) =u(-M,m). (21) 

Finally taking into account that p~ M u(M,m) represents probability intensity transition of the Markov process and 
(21) we may write 

u(-M, m) = a(-M + m - 1) - a(-M + to), (22) 

where {a(— n)}, n = 0, 1, 2, is a sequence of positive numbers such that 

a(-n) > a(-n+ 1) (23) 

Proceeding in the similar way as in [19] we obtain the solution of the Kolmogorov equations for the Markov process 
with local characteristics given by (20) and (21) in the form 

M-l 

pM{t) = P - M + P — V- 1 cxp(tW-M,+i), (24) 
p »=o 

if dist p {Bf, Bf) = p - M + m then 

M-l 

P^\t)^p- M + P~ (m+l) cxp(tW- M ,m +i+ i) -p- m exp(W_M, m ) (25) 

P i=0 

where 

M-l 

W_mj = - J! fc ) - U (~ M > k + l))p~ M+k - (26) 

k=j 

Now we define the transition probabilities of the Markov process on Z p in the following way. 
For any x E Z p let 

P t (x; B(a, M)) =P(X M (t) e S(x, M) | X M (0) G S(o,Af)) , (27) 

By arguments similar to those in [19] we may prove that there exists continuous time Markov stochastic process 
{X(t), t > 0} with state space Z p and transition probabilities Pt(x; B) given by (27). 

Observe that (27) give us direct connection between Markov processes {Xm(<), t > 0} on p~ adic balls with radius 
p~ M , and Markov process {X(t), t > 0} on Z p . 

III. DYNAMICS AS A THERMAL HOPPING IN p ADIC SPACE. 

In this section we shall apply the mathematical considerations from Section 2 to the physical system with the 
hierarchy of states which can be linked with Cayley tree structure. By studying the dynamics of such the systems 
temperature-dependent power law decay and the Kohlrausch law are derived. Transitions between states are thermally 
activated. The height of the energy barriers, k = 1, 2, . . . , which the system overcomes can be ordered in the 
increasing sequence Ai < A2 < • • • < A^ .... The time evolution of the system is described by a random walk on the 
space of states. The simplest model of such dynamics is the model proposed by Ogielski and Stein [9]. They consider 
a regular Cayley tree with M levels and fixed branching ratio p. The total number of leaves, points on the bottom 
of the tree, is n — p M . The natural ultrametric distance d(k, I) between leaves k and I , is defined as equal to the 
height m, m = 0, 1, . . . , M of their closest common ancestor. Now, identifying the states x and y separated by the 
energy barrier A m with leaves k and I, the probability of moving from state x to state y may be defined as equal to 
transition probability from leaf k to leaf I separated by the ultrametric distance to. Thus dynamics in the space of 
states separated by the energy barrier may be studied in terms of appropriate Markov process involving the end points 
of Cayley tree, as a space of states. It is a nontrivial observation that the probability transition intensities between 
states are depending on their ultrametric distance. Due to hierarchical structure of the state space a probability 
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intensities matrix of the process has Parisi matrix structure. Parisi matrix has regular form, and for illustration we 
present the case p = 2. 



ei e 
Ei 



Ei 

ei e 



Ei 



E9 



Ei 

£o e l 



(28) 



where £7, is the matrix with all elements equal to ej. 

One can observe that end points of a regular Cayley tree with M levels and fixed branching ratio p may be represented 
as a set of disconnected balls {Bq 1 , . . . , B^l, _ 1 } with radius p~ M covering Z p . 

Let us consider now a special case of Markov process on p— adic integers. We assume that transition probability 
intensities of the process {X(t), t > 0} depend on p— adic distance only. For this process we have a corresponding 
Markov chain {X M (t), t > 0} with the set of disconnected balls {B^f , . . . , B M Ml } with the radius p M covering 

Z p , as a state space. If we enumerate these balls in such a way that dist p (B Q M , Bf 1 ) increase with i then Q = [qij], 
(0 < j < p M — 1, < i < p M — 1), the matrix of transition probability intensities of the process t > 0} has 

Parisi matrix form. By appropriate choice of u(—M, fc), we obtain process studied in [9]. 
Let for k = 1, 2, M- 1 



p- M u(-M, k) 



q {M) 

Hi] ! 



(29) 



where gjj^ is probability transition intensity of a jump from the ball £? 2 M to the ball Bj 1 , separated by p— adic 



dists 



mce p 



-(M-fe) 



From (24) we have that process which at time starts from a ball Bq 1 will be found at this ball at time t with 
probability 



P t «, B™)=p- M + P —Y. p- i exp(-tW_ M ,i + i)} 



M-1 



i=l 



which together with (26) and (29) gives 



P t (B M ,B M )=p- JM + 



p-1 



M-1 



M 



J^P 4 cxp(-t ^ (e k -e k+1 )p k ) 

k=i+l 



i=0 



For a special case p = 2 equation (31) has the same form as corresponding equation (6) from 



1 M-1 M 

P t (B™, O = 2- M + - 2- < (exp(-t[2a i+1 + ^ a fc ]) 



(30) 



(31) 



(32) 



i=0 



fe=i+2 



where at = 2 fc 1 represents the probability intensity of a jump an ultrametric distance k from a starting sit at 
Cayley tree, while in our case a k represents the probability transition intensity of a jump to any ball of radius 2~ M 
at 2-adic distance 2~( M - fe ). 

Finally we contract the ball radiuses to zero and performing procedure analogical to [9] and specifying the form of 
energy barriers and probability intensities of crossing a barriers we are able to compare different scenario. The most 
simple case is a sequence of barriers linearly growing with p— adic distance = Afc for some positive constant A 
and afc = e~ Ak ' T and for fixed temperature T. In this case (31) gives us fore large time t a temperature-dependent 
power law. 



lim P t {B^,B^)^r Tln2 ' A . 

M— »oo 



(33) 



For a sequence of energy barriers which grows in slower way, i.e Afc — A lnfc, for some positive constant A and 
afc = e~ A lnk / T and for fixed temperature T, the probability P t (B^f , Bq 1 ) fulfils, for large t, the Kohlrausch law 



lim P t «,i? M )~expH T / A ). 

M^oo 



(34) 
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IV. CONCLUDING REMARKS 



In our paper use random walk framework and we show that p— adic analysis is very natural tool to describe the 
relaxation process in a glassy systems, p— adic space has a well defined ultrametric topology which we employ in this 
paper. We do not use the random walk along the Cayley tree [9] because one can measure the ultrametric distance 
between the physical states directly on the bottom of the tree. 

For the description of dynamics of the hierarchical systems besides of the random walk on p— adic numbers one can 
alternatively use the process of jump diffusion on p— adic numbers which description employs the pseudodifferential 
operators [26] in p— adic space. In these both approaches one obtains the equivalent physical results [22-24]. 

It seems that the p— adic analysis is a good tool not only for the description of the dynamics of the glassy systems, 
but also for another hierarchical processes like the evolution of fractals [27], the avalanches [28], protein folding [29] 
ect.. It is interesting to notice that the p— adic space inherently includes the natural hierarchy: p— adic balls can be 
covered by the smaller ones disjoint balls. The hierarchy of these nested balls correspond to the hierarchy of the scales 
of the configuration rearrangements, as it is seen from the scaling theory for the growing domains and droplets [5-7]. 
The droplets may be broken into smaller ones when temperature decrease. 

It is also interesting to notice that by the mapping of any Cayley tree onto its bottom states one can obtain the 
integer p— adic space and further changing the fractional part of the main ancestors we obtain another integer p— adic 
space, p— adic space consists of infinite number of these integer p— adic spaces. In conclusion one can comment that 
the memory effects in spin glassy systems it seems to be are well described by adic topology. 
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